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Several useful thermodynamic relations are derived for metal-insulator transitions, as general- 
izations of the Clausius-Clapeyron and Eherenfest theorems. These relations hold in any spatial 
dimensions and at any temperatures. First, they relate several thermodynamic quantities to 
the slope of the metal-insulator phase boundary drawn in the plane of the chemical potential 
and the Coulomb interaction in the phase diagram of the Hubbard model. The relations impose 
constraints on the critical properties of the Mott transition. These thermodynamic relations are 
indeed confirmed to be satisfied in the cases of the one- and two-dimensional Hubbard models. 
One of these relations yields that at the continuous Mott transition with a diverging charge 
compressibility, the doublon susceptibility also diverges. The constraints on the shapes of the 
phase boundary containing a first-order metal-insulator transition at finite temperatures are 
clarified based on the thermodynamic relations. For example, the first-order phase boundary is 
parallel to the temperature axis asymptotically in the zero temperature limit. The applicability 
of the thermodynamic relations are not restricted only to the metal-insulator transition of the 
Hubbard model, but also hold in correlated systems with any types of phases in general. We 
demonstrate such examples in an extended Hubbard model with intersite Coulomb repulsion 
containing the charge order phase. 

KEYWORDS: charge compressibility, doublon susceptibility, Mott transition, thermodynamic relation 



1. Introduction 

The interplay between the kinetic energy and the 
correlation effect for Fermion systems has been one of 
the most important subjects in the condensed mat- 
ter physics. The interplay of quantum fluctuations and 
many-body effects can induce highly non-trivial phenom- 
ena. The metal-insulator transition driven by the elec- 
tron correlation, which is called the Mott transition, 1 
offers such a prototype. 2 

In the Mott transition, there are two control param- 
eters: one is the chemical potential and the other is 
the bandwidth defined by the relative strength of the 
Coulomb interaction to the hopping integral. Controls by 
these two parameters are realized in a lot of examples in 
real materials as transition-metal compounds 2 including 
cuprates which exhibits the high-T c superconductivity, 
organic materials 3 and 3 He systems. 4 

Recently, critical exponents of the Mott transition 
through the bandwidth-control route in the phase dia- 
gram of temperature, pressure and magnetic field have 
been reported in (Vi-^Cr^Oa. 5 The critical point was 
also identified in k-(ET) 2 Cu[N(CN) 2 ]C1. 6 Although the 
measurements in (Vi-^Cr^Oa support the Ising uni- 
versality 5 predicted by G. Kotliar, et al. 7 ~ 9 within the 
dynamical-mean-field-theory in the infinite spatial di- 
mension, 10 the measurements in k-(ET) 2 Cu[N(CN) 2 ]C1 
suggest a new aspect for the Mott transition. 6 These cir- 
cumstances indicate the importance of the theoretical 
study to clarify the basic properties of the Mott transi- 
tion and the mutual relationship among the basic phys- 
ical quantities near the Mott transition in finite spatial 
dimensions. 

Theoretical studies of the Hubbard model by fully tak- 
ing into account the spatial fluctuations in finite dimen- 
sions have been done by numerical calculations and an- 



alytical methods. For the filling-control Mott transition, 
the continuous character with the diverging charge com- 
pressibility was clarified at zero temperature in the one- 
dimensional Hubbard model by the Bethe-Ansatz calcu- 
lation 11 and in the two-dimensional Hubbard model by 
the quantum Monte Cairo (QMC) calculations. 12, 13 

For the bandwidth-control Mott transition, the 
ground-state phase diagrams of the half-filled Hubbard 
models on the square lattice 14 and the anisotropic trian- 
gular lattice 15 have been determined by the path-integral 
renormalization group (PIRG) method. 16 ' 17 It has been 
suggested that the order of the Mott transitions changes 
from the first-order to continuous ones as the effect of 
magnetic frustration becomes large. 

Recently, a new numerical algorithm, the grand- 
canonical path-integral renormalization group (GPIRG) 
method has been developed, 18 which enables us to study 
the filling- and bandwidth-control Mott transitions in a 
unified framework. The GPIRG method has first made 
it possible to determine the ground-state phase diagram 
in the plane of the chemical potential and the Coulomb 
interaction of the square-lattice Hubbard model. 18 It has 
been found that the V-shaped Mott insulator phase ap- 
pears with the first-order bandwidth-control Mott tran- 
sition at the corner of the V shape and the contin- 
uous filling-control Mott transition at the edges. The 
V-shaped structure has been shown to be consistent 
with the contrast found in the difference between the 
first-order bandwidth- and the continuous filling-control 
routes. This consistency is derived from the generalized 
Clausius-Clapeyron and Ehrenfest relations between the 
slope of the metal-insulator boundary and physical quan- 
tities. 18 

In this paper, we further study the thermodynamic 
relations. We show that a couple of general relations 
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between physical quantities hold at each case of the 
first-order and continuous Mott transitions. These re- 
lations are again derived from the Clausius-Clapeyron 
and Ehrcnfcst relations and hold at any spatial dimen- 
sions and at any temperatures. We derive the thermody- 
namic relations between the charge compressibility and 
the doublon susceptibility for the continuous Mott transi- 
tion, where the doublon susceptibility is shown to diverge 
when the charge compressibility diverges. The relations 
impose constraints on the shapes of phase diagrams. 

The organization of this paper is as follows: In §2, 
the basic physical quantities for filling- and bandwidth- 
control routes are defined. The thermodynamic relations 
between the slope of the phase boundary and the physi- 
cal quantities are derived for each case of the first-order 
and continuous Mott transitions. In §3, it is shown that 
they are actually satisfied in analytic forms obtained in 
the one- and two-dimensional Hubbard models. In §4, 
the thermodynamic relations are extended to finite tem- 
peratures. The constraints on the shapes of the phase 
diagram for the finite-temperature first-order transition 
are classified by using the thermodynamic relations. In 
§5, an application to other systems is demonstrated by 
taking an example of an extended Hubbard model with 
intersite Coulomb repulsion. We summarize the paper in 
§5- 

2. Derivation of thermodynamic relations 

2. 1 Basic physical quantities for filling- and bandwidth- 
control routes 
The Hubbard model which we consider is 

H = - *y ( c L c i<? + c )a c i<y) - M^ n ^ 



+u 



N 
i=l 



(1) 



where da (cJ CT ) is the annihilation (creation) operator on 
the i-th site with spin a and n icr = c\ a c ia in the N- 
lattice system. Here, Uj is the transfer integral and in 
this paper, the nearest-neighbor transfer is taken as the 
energy unit. In cq. (1) \i is the chemical potential and U 
is the Coulomb interaction. 

The electron filling n and the double occupancy D are 
defined by the first derivative of the Hamiltonian as 



n 



D 



N d/j, ~ t\l 2^"*>' 



(2) 



i— la 



i d(H) 

N dU 



2—1 



(3) 



where (• • • ) represents the ground-state expectation 
value or the thermal average at finite temperature. 

The charge compressibility Xc and the doublon sus- 
ceptibility xd arc defined, respectively, by the second 



derivative of the Hamiltonian as 
1 d 2 (H) 

Xc ~ n dfi 2 
i d 2 (H) 



Xd 



N dU 2 



dD 

dU 



(4) 
(5) 



2.2 The case where the phase boundary U(n) has finite 

slope SU/S^i 
2.2.1 First- order-transition case 

Although the following derivation of eq. (8) has been 
given in ref., 18 we show here it for the self-contained de- 
scription. Let us consider the ground state, although the 
results described below hold at finite temperatures, as 
will be mentioned in §4. The ground-state energy is ex- 
panded by \i and U: 



E(n + §n,U + SU) 
= E(p,U)- 



dU 



su 



+0(6n 2 ,5U 2 ). 



(6) 



In the fi-U phase diagram, along the metal- insulator- 
transition boundary, the energies of the metal and in- 
sulator phases should be the same because of the coexis- 
tence: 

E^U) - E M (ji,U) 

= Ei{fi + S^,U + SU) - E M {n + S^,U + SU) 



= 0, 



(7) 



where subscripts I(M) represent that the expectation 
value is taken in the insulating (metallic) ground state. 
If the first-order metal-insulator transition occurs, the 
slope of the transition line in the ji-U phase diagram is 
determined by the ratio of the jump of filling and double 
occupancy. Namely, from eq. (6) and eq. (7) the following 
equation is derived: 

SU _ ni - n M 
5n~ Di-Dm' 

2.2.2 Continuous-transition case 

Let us consider the phase diagram in the fj,-U plane 
in the ground state, T = 0. If the continuous metal- 
insulator transition occurs and the slope of the transition 
line is finite in the /i-U plane, the double occupancy is 
expanded by \i and U, as follows: 



D(n + 5ii,U + 5U) = D{pi,U) 

SU 



dD\ . (dD 



dU 



+ 0(S^ 2 ,SU 2 ).(9) 



Along the metal-insulator-transition boundary, the dou- 
ble occupancies at the metallic and insulating phases 
have the same value at (fi, U) and (fi + Sfi, U + SU): 

D M {n + S^,U + SU) - D 7 (/i + Sfi, U + SU) 

= D M (ii,U)-D I ( i i,U), (10) 

where the subscripts M(I) represent that the average 
(• • • ) in eq. (3) is calculated in the metallic (insulating) 
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state. By substituting eq. (9) to eq. (10), the slope of the 
metal-insulator-transition boundary is expressed as 



SU 



dD \ 

auJ, 



M 



_ (9D) 



(od\ 




Ju 


M 



(11) 

Xd\ m - Xd^ 

where {dD/d^)u\i = and eq. (5) are used. 

On the other hand, we can derive another expression 
of SU/S/i starting from the expansion of filling n. 18 In 
the parallel discussion as above, the slope of the metal- 
insulator transition boundary is derived as 

SU 
Sfj, 



Xc 



\dU) , 



(12) 



If (dn/dfi)u and (dD/dU)^ are continuous at (//, U), the 
following relation holds: 



On 
dU 



We then obtain 
dn 



1__5_ f d(H) \ 
NdU V dfi J u 

1 g (d{H)\ 
Ndfi\ dU J ^ 



3D 



dU 



M 



3D 



Xc(x 



DIM 



XD\I) 



(13) 



(14) 



M 



By applying this equation to the thermodynamic rela- 
tions of eq. (11) and eq. (12), the following Theorem 1 
is derived: When the metal-insulator transition is a con- 
tinuous one with a finite and nonzero slope in the ji-U 
phase diagram, the divergence of the charge compress- 
ibility occurs simultaneously with the divergence of the 
doublon compressibility. If the slope of the phase bound- 
ary is zero, SU/S/j, = 0, \c —> 00 Xd\m —> °o still 
holds, whereas the converse is not necessarily true. 

It should be noted that Theorem 1 derived above 
holds irrespective of the spatial dimensionality. The con- 
ditions in which Theorem 1 holds are as follows: First, 
there exists a finite length of the continuous metal- 
insulator-transition boundary U (/i) in the fi-U plane, 
where SU(fj,)/Sfj, exists. Second, the slope of the metal- 
insulator boundary is finite, i.e., \SU/Sfi\ < oo. If the 
metal-insulator transition takes place as the continuous 
one with the U-shaped structure as shown in Fig. 1(a), 
the thermodynamic relation holds at any point of the 
transition line. 

2.3 The case where the slope SU/Sfi is not well defined 
When the Mott-insulator phase has the V-shaped 
structure in the fjrU phase diagram as shown in Fig. 1(b), 
the above discussion is not applied at the corner of the V- 
shaped boundary, which is the bandwidth-control Mott- 
transition point, since 5U/Sfx is not well defined at the 
point. The V-shaped Mott-insulator phase in the \i-U 
plane seems to appear actually in the Hubbard model 
on the square lattice for t'/t = —0.2 with t and t' be- 
ing the nearest- and next-nearest-neighbor transfers. 18 



(a) 
i insulator ; 




U 



metal 



(b) 
. insulator/ 



U, 




metal 



(c) 
insulator. 




metal 



Fig. 1. Classification of the shapes of the phase diagrams in the 
plane of the chemical potential /i and the Coulomb interaction U: 
(a) U shape, (b) V shape, and (c) T shape. The solid lines repre- 
sent the continuous metal-insulator transition. In (b) and (c), the 
open circle located at the corner indicates the bandwidth-control 
Mott transition point to which the thermodynamic relation can- 
not be applied (see text). 



In this case, the first-order bandwidth-control Mott tran- 
sition at the corner is compatible with the continuous 
filling-control Mott transition at the edges. This offers 
an example that the corner of the V shape can have the 
different character of the transition from the edges. 

In the one-dimensional Hubbard model without 
the next-nearest-neighbor hopping, the metal-insulator 
phase boundary has a shape of fi ~ VtUexp(~2-irt /U) 
for the sm&W-U/t regime. 19 In the Hubbard model on 
the square lattice without the next-nearest-neighbor hop- 
ping, the metal-insulator phase boundary is expected to 
have a shape of fi ~ texp(—2ny/t/U) for the small- U/t 
regime. 20 Both systems with the perfect nesting have the 
essential-singular form of the phase boundary, which is 
classified into the T-shaped case as shown in Fig. 1(c). 
This is categorized to the same class as the V-shaped 
case. Namely, the thermodynamic relation can be applied 
except for the corner of the T-shaped boundary. 

3. Analytic form of doublon susceptibility 

In this section, we derive some useful analytic forms of 
dn/dU, dD/d/j, and \d = —dD/dU near the continuous 
Mott transition, starting from the analytic forms of Xc 
realized in the one- and two-dimensional Hubbard mod- 
els, and examine the thermodynamic relations, eqs. (11), 
(12) and (13) derived in §2. In this section, we approach 
the transition point from the metallic side. 

3.1 Derivation in one- and two-dimensional Hubbard 
models 

Let us consider the vicinity of the metal-insulator 
transition in the metallic side of the fi-U phase di- 
agram. In the one-dimensional Hubbard model with- 
out next-nearest-neighbor hopping, 11 and in the two- 
dimensional Hubbard models with 13 ' 18 and without 12 
the next-nearest- neighbor hoppings, the filling depen- 
dence of the chemical potential is proposed to be scaled 
as 

H = fi c - -^S 2 . (15) 

with 8 = 1 — n and a constant a. In the one-dimensional 
system, eq. (15) was derived by the Bethe-Ansatz anal- 
ysis 11 and in the two-dimensional system, eq. (15) was 
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Fig. 2. Two ways of approaching the metal-insulator transition 
point denoted by the black circle at (ii c ,U c ) in the fi-U plane: 
One is U — > U c — at fi = fi c and the other is fi — > fi c — 
at U = U c - These are the routes from the metallic side. The 
shaded area is the insulating phase and the solid line is the phase 
boundary. 



identified by the QMC calculations. 12, 13 Note that here 
we follow the notation of the chemical potential in 
refs. 12,13 to facilitate the comparison. From eq. (15), the 
doping rate is expressed by the chemical potential as 



5 = ^2a(fi c — /a). 



(16) 



We need to consider the electron doping separately when 
the next-nearest neighbor hopping becomes nonzero 
causing the electron-hole asymmetry. Here, we consider 
the hole-doped case, although the following discussion 
can easily be extended to the electron-doped case. Then, 
the filling dependence of the charge compressibility is 
given by 

Xc = z - T - 



1-71 



By substituting eq. (16) to eq. (17), the chemical- 
potential dependence of the charge compressibility is 
given by 



Xc 



|(Mc 



m) 



-1/2 



(18) 



To facilitate the discussion below, the most divergent 
terms of dD/d^i, dn/dU and xd = —dD/dU near the 
metal-insulator transition are characterized by the expo- 
nents p, q and r defined below. 

We consider the metal-insulator phase boundary illus- 
trated in Fig. 2. Two routes approaching the transition 
point at (/i c , U c ) indicated by the two arrows are now 
the subject of the present study: One is Z7 — > U c — at 
[i = fi c and the other is fi — ► /i c — at U = U c . 

At U = U c , the relation between n and fj, was already 
given in eq. (16), while the double occupancy as a func- 
tion of fj, near \i c may be similarly assumed as 

c( Mc - M ) 1 -p 



D 



for fi < He, 



Dhf r rr- .,. 

D hf for Mc < M <0. 1 



Here Dhf is the double occupancy at half filling (in the 



one-dimensional Hubbard model, Dh / is analytically cal- 
culated in ref. 21 ). 

The most dominant term for the filling as a function 
of U near U c is assumed in the following form: 



-a(U c - U) 1 - 11 + 1 



1 



for U < U c , 
for U c < U. 



(20) 



At [i — fj, c , the double occupancy as a function of U 
near U c may be assumed in the following scaling form: 

_ f d(U c - Uf- r + D hf (U+) for U <U C 



D = 



D hf {U) for U C <U 

In the metallic phase, the following equations are de- 
rived from eq. (19) and eq. (20), respectively: 

(<9m);7 



(21) 



-c{l-p)(n c -n)-v, (22) 



dn 
dU 



M 



M 



= a(l-q){U c -U)- q . 



(23) 



From eq. (21) the doublon susceptibility is written as 
Xd\ m = d(l - r){U c - U)- r for U < U c , (24) 
9D hf (U) 



dU 



for U c < U. 



(25) 



In the following, the coefficients a, c, and d, and the 
critical exponents p, q and r are derived analytically. Let 
us first consider eq. (19). From eq. (15), the ground-state 
energy is obtained by integrating fi over n: 

E(S)=E hf -fi c 5+-^8 3 . (26) 



6a 



From this equation, the double occupancy is obtained as 
„ dE fdn\ . 1 fda\ r3 

D = m = D "-{m)J-**[m)f- (27) 

Now we consider the vicinity of half filling in the metallic 
phase and expand D up to the first-order in term of 8. 
By substituting eq. (16) to eq. (27), D is expressed by fi: 



(28) 



for the most dominant term. By differentiating D by fi, 
we obtain 



dD 



a ( dfi 
2 { dU 



(29) 



for the most dominant term. Comparing eq. (29) with 
eq. (22), we obtain the relation between the coefficients 
c in eq. (19) and a: 



dU 



(30) 



and the critical exponent is obtained as p = 1/2. 

Second, let us consider eq. (20). Starting from eq. (16), 
we obtain 

d6_ 

dU , 



+ 



!(*-,)-*(*) pi) 
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Fig. 3. Two ways of approaching the metal-insulator transition 
point denoted by the black circle at (/i C; Uc) in the fi-U plane: 
One is U — > U c — at fi = fi c and the other is fi — » fjb c — 
at U = U c - These are the routes from the metallic side. The 
shaded area is the insulating phase and the solid line is the phase 
boundary. 



in the metallic side around (/i c , U c ). Here, we consider the 
point (/i, U) close to the phase boundary as in Fig. 3. If 
we define the slope of the phase boundary near (/i, U) as 
(d/j,/dU) c , the following relation is obviously satisfied: 

d\x ' 
dU , 



{fl>c ~ (J>)\u= Uc 



(32) 



In eq. (31), in the limit of /i — > fi c — 0, the second term 
is dominant and by using the relation, eq. (32), the di- 
verging part is written as 

1/2 

{U c -U)- 1/2 . (33) 



On 
dU 



d5_ 

dU 



dU 



Then, the coefficient a in eq. (20) and eq. (23) is ex- 
pressed by a and the slope of the metal-insulator bound- 
ary: 

a = ^{-^j'\ (34) 

and the critical exponent is obtained as q = 1/2. 

Third, let us consider eq. (21). By the derivative of 
both sides of eq. (28) by U, the doublon susceptibility is 
written as 

)2, 



dD 

dU 



dD 



dU 



1/2 



dU 

dfi 
dU 



'2a(fi c - fi) 



w) (/ic 



1/2 



(Mc-M)- 1/2 - 



(35) 



Here (d 2 fi/dU 2 ) c is the curvature of the phase boundary. 



In the limit of /i 



the most dominant term is the 



third term of the right-hand side of eq. (35), and by using 



eq. (32), we obtain 



3/2 



Then, the coefficient d in eq. (21) and eq. (24) is ex- 
pressed by a and the slope of the metal-insulator bound- 
ary: 

3/2 

d=V2a[-^\ , (37) 



and the critical exponent is obtained as r = 1/2. 

By eqs. (30), (34) and (37) we have derived the dom- 
inant terms of (dD/d/j.) | u=u c , (9n/9C/)| M = Mc and \d = 
— (dD / dU)\jj,=n c near the Mott transition, starting from 
the singularity of \c- Next, let us discuss the relationship 
between the coefficients of these terms and the slope of 
the metal-insulator boundary, and test Theorem 1. 

3. 2 Examination of the thermodynamic relations 

From eqs. (30), (34) and (37), the following relation is 
derived: 

mr 



2a 



a f a \' 
d = \c) 



(38) 



So far, it is shown that the coefficients a, c and d are 
expressed by a and (dfi/dU) c as in eqs. (30), (34) and 
(37), and the relation of eq. (38) is satisfied. To examine 
the thermodynamic relations of eq. (11) and eq. (12), 
let us focus on eq. (13): In the metallic phase, dn/dU 
and dD/d/j, are continuous and hence eq. (13) holds. In 
the vicinity of the metal- insulator boundary, dn/dU and 
dD/dfi have the form of eq. (22) and eq. (23) and the 
equation holds: 

a(U c - U)- 1 / 2 = c{n - Me)" 172 . (39) 

By using the relation of eq. (32), the following limit is 
taken in eq. (39): 



lim c 



lim a(U c - Uy 1/2 



(40) 



We see that this equation is satisfied by eq. (38). This 
implies that eq. (13) holds in the limit to the metal- 
insulator boundary from the metallic phase in the fi-U 
plane as it should be 

m)- 1/2 - 



lim a(U c - U)- 1/2 = lim c(^ c 
By using eq. (41), eq. (38) is rewritten as 



(41) 



dU 



lim. 



2a(fi c - fi)- 1 ' 2 



lim 



^{Uo-U)- 1 ' 2 ' 



lim,,..,- c(u c - u) 1 / 2 



lim 



d(U c -U)~V 2 ' 



This demonstrates how the thermodynamic relations of 
eq. (11) and eq. (12) with a simultaneous diverging of 
the charge compressibility and the doublon susceptibil- 
ity hold. Then, the validity of the thermodynamic re- 
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Fig. 4. Ground-state phase diagram of the one-dimensional Hub- 
bard model for t = 1. The metal-insulator phase boundary is 
obtained by the Bethe-Ansatz solutions. 19 



lations is confirmed analytically in the case with Xc — 
\/a/2(/i c — Z- 4 ) -1 ^ 2 which is known to be realized in the 
one- and two-dimensional Hubbard models. 

3.3 Bethe-Ansatz solutions in one- dimensional Hubbard 
model 

In the one-dimensional system, the exact solutions of 
the Hubbard model are available by the Bethe-Ansatz 
equations. 11, 19:21 The energy and the double occupancy 
at zero temperature for the given filling and Coulomb in- 
teraction were calculated in ref., 21 and the coefficient of 
the charge compressibility a in eq. (15) was obtained 
in ref. 11 The chemical potential at which the metal- 
insulator transition occurs for U is given by 



2 r 



1 



dw- 



Jx(w) 



(43) 



where upper(lower) sign denotes the elect ron( hole) - 
doped regime and J\{w) is the first Bessel function. 19 
The ground-state phase diagram of the one-dimensional 
Hubbard model is shown in Fig. 4. By using the exact 
expressions of D in eqs. (2.14) and (2.15) in ref. 21 and 
n in eqs. (6a) and (6b) in ref., 11 we plot the critical be- 
haviors of these quantities as functions of /x and U near 
the Mott transition. As a typical example, we show the 
critical behaviors of n and D around the Mott transition 
point, (/i c , U c ) = (-3.277, 10.0) (see Fig. 4). 

The chemical-potential dependence of the double oc- 
cupancy for U — U c is shown in Fig. 5. We also plot 
D = D h f + \Z2a(—d/j,/dU) c (fi c — /i)~ 1//2 obtained from 
eqs. (28) and (30) by the dashed line. We see that dD/dfi 
diverges at /j = /i c , whose critical behavior is described 
by eq. (29). 

The interaction dependence of the filling for fi = 
[i c is shown in Fig. 6. We also plot n = 1 — 
V2a{-d^/dU)\ /2 {U c - U) 1 ' 2 obtained from eqs. (20) 



-0.21-5 



-0.2S-0 



-0.22-5 




Fig. 5. Double occupancy D vs. chemical potential fj, obtained 
by the Bethe-Ansatz solutions of the one-dimensional Hubbard 
model (solid line) for t = 1 at U = U c = 10. Dashed line repre- 
sents D = D hf + V2a{-d^/dU) c {n c - fJ.)~ 1/2 - 



t=i.o |^=|V=-3.2737 



1.0 3- 



c . 9 



0.93 




10 . 01 



Fig. 6. Filling n vs. Coulomb interaction U obtained by the 
Bethe-Ansatz solutions of the one-dimensional Hubbard model 
for t = 1 at fj, = fM c = —3.2737. Dashed line represents 
n = 1 - V2^(-dfJ./dU)l /2 (U c - U) 1 / 2 . 



and (34) by the dashed line. We see that dn/dU di- 
verges at U = U c , whose critical behavior is described 
by eq. (33). 

The interaction dependence of the double occupancy 
for pk = fi c is shown in Fig. 7. We also plot D = Dhf + 
V2a(-dn/dU)l /2 (U c - C/) 1 / 2 obtained from eqs. (21) 
and (37) by the dashed line. We see that the doublon 
susceptibility diverges at U — U c and its critical behav- 
ior is described by eq. (36). 

In the one-dimensional Hubbard model, the limits of 
a in eq. (15) were shown as a — > for U — * and 



a -> 1/(2tt 2 ) for U - 
C/^Oand {dn/dU) c - 
and d are as follows: a 
and a — > l/(v2~7r), c - 
U -> oo. 



34 



Since (dfj,/dU) c 



for 



1/2 for U —> oo, the limits of a ,c 
> 0, c -> 0, and d -> for U -> 0, 
• 1/(2tt) and d -> 1/(2 V2tt) for 



General expressions in d-dimensional systems with 
dynamical exponents z 
We have discussed the analytic forms of Xc\u=u c , 
(dn/dU)]^^, (dD/dfi)\ u=Uc and Xd\h=h b , which arc 
realized in the one- and two-dimensional Hubbard mod- 
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Fig. 7. Double occupancy D vs. Coulomb interaction U obtained 
by the Bethe-Ansatz solutions of the one-dimensional Hubbard 
model for t = 1 at fj, = /Lt c = —3.2737. Dashed line represents 
D = D hf + V2^{-dix/dU)l /2 {U c - t/) 1 / 2 . 



els. Here we present the general expressions of these 
quantities in d-dimensional systems. 

Based on the hyperscaling analysis, the critical be- 
havior of physical quantities near the continuous filling- 
control Mott transition was predicted in ref. 22 For exam- 
ple, the charge compressibility has the form of 



Xc 



(M C -Mr (d+z) - 2 , 



(44) 



where v is the correlation-length exponent which satisfies 
vz = 1 with z being the dynamical exponent, and d is 
the spatial dimension. It is known that for the usual one- 
dimensional Hubbard model, where d = 1, the dynamical 
exponent is given by z = 2, then the charge compressibil- 
ity has the form of Xc ~ (a*c — f) 1 ^ 2 as in eq. (18). For 
the two-dimensional system, d — 2, the dynamical expo- 
nent z = 4 was identified by numerical calculations, 23,24 
and then the same critical behavior for Xc as the d = 1 
case appears. For the three-dimensional system, in the 
case of z = 4, the charge compressibility has the form of 
Xc ~ (Mc — m)~ 1//4 by eq. (44). In the following, we show 
the analytic forms of the physical quantities defined by 
the second derivative of the Hamiltonian in general form, 
which are derived from eq. (44): 

v{d+z)-l 



XD 



dn 
dU 



d(J>/i-i 



dD 

dU 



v(d+z)-2 



dU 

x(U c -U) 

a ( dfi 
2 \dU 



x( Mc - M r^>- 2 , 

2 \dU 



x{U c ~U) u{d+z) - 2 . 

The coefficients of each term are expressed by the coeffi- 
cients of Xc and the slope of the metal-insulator bound- 
ary. The exponents are the same as that of the charge 



compressibility, i.e., v(d + z) — 2. Note that the relation 
eq. (38) still holds in this case, and this implies that 
eq. (38) does not depend on the exponent in the diver- 
gence of the charge compressibility. This is ascribed to 
the fact that thermodynamic relations of eq. (11) and 
cq. (12) hold in any spatial dimensions. 

4. Finite temperature 

4-1 Generalization of thermodynamic relations to finite 
temperature 

In §3, the thermodynamic relations are examined at 
zero temperature. It should be noted, however, that the 
thermodynamic relations of eqs. (8), (11) and (13) hold 
even at finite temperature, if (...) is taken as the thermal 
average: 



(A) = tr ( G ~ H/T ^ 



(45) 



tr ( e -ff/T) ' 

where A is an operator for a physical quantity and T is 
temperature. Hence, when temperature is fixed, in the 
plane of /x and U with the finite length of the metal- 
insulator boundary, eqs. (8), (11) and (13) hold. 

4-2 Constraints on the shape of the phase diagram for 
finite-temperature first-order metal-insulator tran- 
sition 

We also discuss the phase diagram with the temper- 
ature dependence. Starting from the expansion of the 
Gibbs free energy G(T, U) for electrons along the first- 
order mctal-insulator-transition line, we obtain 



SU 



Si — Sm 



(46) 



Dj - D M ' 

where the entropy of electrons is defined by S = 
— (dG/dT)u and D is the double occupancy in the ther- 
mal average. In experiments, when the pressure p is ap- 
plied to the system, the bandwidth of the electronic sys- 
tem, U/t changes. Hence, in eq. (46), U corresponds to 
p and the double occupancy D is related to the conju- 
gate variable, the volume v. The similarity between the 
metal-insulator transition and the liquid-gas transition 
was pointed out in ref. 25 

We discuss the consequence of eq. (46). Let us con- 
sider the case that the first-order metal-insulator tran- 
sition takes place in the U-T phase diagram and the 
insulator phase has the ordering such as the magnetic 
one. In the metallic phase, at small but nonzero particle- 
hole excitations induce the entropy stemming from the 
finite density of states at the Fermi level, while in the 
insulator phase the entropy is rather small if the symme- 
try breaking with small degeneracy is realized while the 
particle-hole excitation cannot induce the entropy under 
the finite charge gap. Then, this yields Sj — Sm < and 
by using the fact that Di — Dm < 0, the slope should be 
ST/SU > (see Fig. 8(a)). This is a reasonable result, 
since the metallic phase appears in the high-temperature 
side to govern the — TS term in the free energy. This sit- 
uation can be realized in the three dimensional system, 
since Mermin- Wagner's theorem 26 allows the magnetic 
ordering at finite temperature. 

In the two-dimensional system, the magnetic ordering 
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Fig. 8. Some examples of schematic phase diagrams with first- 
order metal-insulator transition in the plane of temperature T 
and Coulomb interaction U. (a): Paramagnetic metal and an- 
fifcrromagnetic insulator phases in the three dimensional sys- 
tem, (b): Paramagnetic metal phase for T > 0, the anfifcrro- 
magnetic insulator phase at T = and the paramagnetic in- 
sulator phase for T > in the two dimensional system. Open 
circle in (b) represents critical end point, T cr , of the first-order 
transition line, (c): Phase diagram in the plane of double occu- 
pancy D and entropy S corresponding to (b). Gray area rep- 
resents the phase-separated region between metal and insulator 
along the first-order-transition line in (b). Dashed line represents 
each set of (T,U) on the first-order-transition line in (b). (d): 
Schematic phase diagram for the continuous transitions both at 
T = and T = T cr drawn by open circles, while the first-order 
transition occurs for < T < T cr . In this case, a finite slope 
{dT / 8U)\t— >o ^ oo is allowed in contrast to (a) and (b)(see 
text). 

at finite temperature is suppressed by thermal fluctu- 
ation. 26 In this case, the paramagnetic insulator phase 
can be realized for finite-temperature regime above the 
magnetically-ordered insulator phase at T = 0. Although 
it depends on each system how the entropy is released by 
the inter-site interaction in the paramagnetic-insulator 
phase, we refer ref. 27 as an example of the appearance 
of the paramagnetic insulator phase at high-temperature 
side, i.e., 5T/6U < by Sj - S M > (see Fig. 8(b)). 

Let us focus on the slope of the first-order transition 
in the zero-temperature limit; (dT/dU)\T-+o- When the 
first-order metal-insulator transition occurs at T = 0, 
the jump of the double occupancy appears at U = U c , 
i.e., D[ — Dm 0. By the requirement of the thermo- 
dynamic third law, the entropy at T = should be in 
the metal and insulator phases, and hence Si = = Sm 
at T = 0. Therefore, from eq. (46), the following con- 
clusion is derived: In the phase diagram with tempera- 
ture T and control parameter U , the first-order transition 
line should have the infinite slope in the zero-temperature 
limit: (dT / 8U)\t^o = oo. This statement is not re- 
stricted to the case of the metal-insulator transition, but 
can be applied to the first-order transition between any 
types of the phases. The schematic phase diagram in the 
plane of the double occupancy D and entropy S cor- 
responding to Fig. 8(b) is illustrated in Fig. 8(c). The 
gray area represents the phase-separated region between 
the metal and the insulator along the finite-temperature 
first-order transition line in Fig. 8(b). The dashed line 



Author Name 

in Fig. 8(c) indicates each set of (T,U) on the finite- 
temperature first-order transition line in Fig. 8(b). 

Here, the relationship to the known results should be 
mentioned. In the infinite-dimensional Hubbard model 
with the fully frustrated lattice, the T-U phase dia- 
gram was obtained: 28,29 The first-order transition be- 
tween the paramagnetic metal and the paramagnetic in- 
sulator takes place at finite temperature. One might con- 
sider that the first-order-transition line has a finite slope 
as (dT / 8U)\t^o < in the phase diagram, which is 
different from the above statement. However, at T = 
the metal-insulator transition takes place as the contin- 
uous one in that model, to which the above discussion 
cannot be applied. Actually, in this case, the continuous 
transition with Di — Dm = occurs at T — 0. Then, the 
finite slope is possible; \dT / 8U\t^vi < oo. The schematic 
phase diagram in the U-T plane for such a case is illus- 
trated in Fig. 8(d): At T = and T = T cr , the con- 
tinuous metal-insulator transitions take place, while the 
first-order transition occurs for < T < T cr . 

At T = 0, possible shapes of the phase diagram with 
the first-order and continuous quantum phase transitions 
have been classified in ref. 18 In the plane of the chemical 
potential fi and the Coulomb interaction U, the U-shaped 
phase boundary with a finite length of the first-order 
transition can exist, but the V-shaped phase is allowed 
to have the only first-order-transition point at the corner 
and have the continuous transitions at the edges. The T- 
shaped phase as an essential singular form of the Mott 
gap, ~ exp[— at/(U c — U)} is classified into the same class 
as the V-shaped case. 

For the U-shaped case at T = 0, a possible T-U-fi 
phase diagram is illustrated in Fig. 9(a). The surface of 
the first-order transition is formed as drawn by the gray 
area in Fig. 9(a) and at the critical end curve drawn by 
the solid line, which describes T cr , the charge compress- 
ibility diverges. 30 For the V-shaped case at T = 0, the 
possible T-U -[i phase diagram is illustrated in Fig. 9(b). 
In this case, at T = the first-order transition occurs at 
the corner of the V-shaped boundary. At finite tempera- 
tures, one possibility is that the jump of the double occu- 
pancy at T = disappears for infinitesimal temperature. 
Actually one can consider such a function for D(hq, U, T) 
with fio the chemical potential at half filling: An example 
is the functional form analogous to the Fermi distribution 
function, f(e,T) = l/[exp((e - fi)/T) + 1]. At T = the 
jump of f(e,T) at e = fj, appears, while at infinitesimal 
T, the jump disappears. This is the case illustrated in 
Fig. 9(b) where T cr becomes zero temperature, i.e., the 
quantum tricritical point. In this case, D(^lq,U,T = 0) 
as a function of U has a jump at the bandwidth-control 
Mott transition point, U c q, and the jump disappears at 
U c0 for T ^ 0. 

The other possibility is that the jump in D(/iq,U,T) 
remains even for T ^ 0. In this case, the first-order- 
transition line, but not the surface as in Fig. 9(a), has 
a certain length from T = to T = T cr . This seems 
to be unusual, but such a case is not ruled out by the 
thermodynamic argument. 
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surface of 




Fig. 9. Possible schematic phase diagrams in the parameter space 
of temperature T, chemical potential /j and interaction U for (a) 
U-shaped and (b) V-shaped Mott-insulator phases at T = 0. In 
(a), thick(thin) line represents the first-ordcr(continuous) tran- 
sition line. Gray surface represents the surface of the first-order 
transition and its upper edge is the critical temperature, T cr . In 
(a) and (b), the solid circle represents the band- width control 
Mott transition point at U c o, and /irj represents the chemical po- 
tential at half filling at U = U c o ■ The dashed line connected to fio 
is the route of the band- width control. In (b), the first-order tran- 
sition takes place only at the bandwidth-control Mott-transition 
point, i.e., T cr = 0. 



5. Application of thermodynamic relation to 
general models 

In §2, we have derived the thermodynamic relations 
for the metal-insulator transition, and in §3 we have 
examined that physical quantities defined by the sec- 
ond derivative of the Hamiltonian for fi and U actually 
satisfy the thermodynamic relations in the cases of the 
one- and two-dimensional Hubbard models. However, the 
thermodynamic relation between the slope of the phase 
boundary and physical quantities defined by the first and 
second derivative of the Hamiltonian for the first-order 
and continuous transitions, respectively, can be derived 
in general in the other models with any types of phases 
in the same way as in §2. 

As an example, let us consider the phase diagram 
whose control parameters are the on-site Coulomb inter- 
action U and the nearest-neighbor Coulomb interaction 
V. The Hamiltonian is described by the extended form 
as 



H 



(47) 



(id) 



where denotes the nearest-neighbor sites. 

At quarter filling, the existence of two phases, namely, 
metallic and charge-ordering phases has been reported in 
the one-dimensional system. 31 

At the first-order transition, by expanding the ground- 
state energy by U and V in a similar way to the deriva- 
tion of eq. (8), the slope of the phase boundary is ex- 
pressed as 

SU R\ M - RU (48) 



sv 



D 



M 



D 



where 



R 



1 d(H) 
N dV 



(id) 



Here, at the first-order transition, the slope of the phase 
boundary in the U-V plane is expressed by the ratio of 
the jump of the nearest-neighbor correlation and the dou- 



ble occupancy. 

At the continuous transition, the slope of the phase 
boundary is expressed in two ways: One is the expansion 
of the double occupancy and the other is the nearest- 
neighbor correlation. The resultant thermodynamic rela- 
tions are 

Xd\ m 



6V 
SU 



XD\i 




(49) 



Xv\ M ~ Xv\i 

where the susceptibility for the nearest-neighbor correla- 
tion is defined by 

d 2 (H) 



Xv 



dV 2 



(50) 



If (dD/dV)u and (dR/dU) v are continuous at (U,V), 
the following relation holds, which is shown in the same 
way as in the derivation of eq. (13): 

'8R\ , . 

„ v«?)v- <51) 

When we add the axis of the chemical potential to the 
phase diagram in the plane of U and V, the thermody- 
namic relations in the fi-U plane and the \i-~V plane can 
be derived as in a parallel way to §2. 

Here, we show the thermodynamic relations for the 
first-order and continuous transitions in the extended 
Hubbard model as an example. It should be noted that 
the thermodynamic relations can also be applied to other 
systems in a similar way. Hence, we remark that the ther- 
modynamic relations derived in this paper have general- 
ity. The applicability of the relations is not restricted to 
the specific model, but to any types of the models with 
the phase diagram. This is a consequence of the exten- 
sions of the Clausius-Clapeyron and Ehrcnfcst relations 
to microscopic models. 

6. Summary 

We have analytically derived several useful thermo- 
dynamic relations of the charge compressibility and the 
doublon susceptibility to the slope of the metal-insulator 
phase boundary in the plane of the chemical potential /i 
and the Coulomb interaction U. By using this relation, 
we have shown that the charge compressibility divergence 
at the metal-insulator transition causes simultaneous di- 
vergence of the doublon susceptibility when the metal- 
insulator transition takes place as a continuous one with 
a finite length and a finite and nozero slope in the fi-U 
plane. This statement holds at any spatial dimensions 
and at any temperatures. The main results of the ob- 
tained relations are eqs. (11), (12) and (14). 

In one- and two-dimensional Hubbard models, critical 
divergence of the charge compressibility at the Mott tran- 
sition appears generally with the form of Xc = \[%{pc — 
/i) -1 / 2 . Starting from this critically divergent form of \c, 
analytical expressions of {dD/dfi)u—i/ c , (dn/dU)^ =flc 
and Xd|/j=p c 111 the vicinity of the Mott transition are 
derived. All the coefficients of the most divergent terms 
of these quantities are expressed by the coefficient of the 
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charge compressibility, a, and the slope of the metal- 
insulator boundary, (SU/Sfj,) c . All the exponents of the 
divergence of these quantities have the same value as the 
exponent of \c'- namely, —1/2. The relations among the 
coefficients of Xc\u=u c , {dD / 'dfi)u=u c , {dn/dU)^ =lla and 
Xd\h=h c at the metal- insulator-transition point actually 
satisfy the thermodynamic relations derived above, i.e., 



(dny 

\du) 



M 



Xc(xd\m - Xd\i)- 



M 



The above statement holds at any points of the U- 
shaped metal-insulator boundary and at the points of the 
V (T)-shaped boundary except the corner, which is the 
bandwidth-control Mott transition point. At the corner 
of V- and T-shaped phases (5U/Sfi) c is not well defined. 

By applying the above argument to the phase diagram 
in the parameter space of the temperature T and the in- 
teraction U, constraints on the phase boundary are clas- 
sified at finite temperatures. It is shown that the infi- 
nite slope of the first-order-transition line has the infinite 
slope as (dT/dU)\T^o = oo from the requirement of the 
thermodynamic third law and thermodynamic relations. 

These conclusions described above are not restricted 
only to the Hubbard model, but also can be applied to 
other general models. As an example, we demonstrate 
the derivation of the thermodynamic relations in the ex- 
tended Hubbard model with the next-near-neighbor re- 
pulsion. 
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